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^ ■ The two-dimensional hydrogen with a hnear potential in a magnetic field is solved by two dif- 
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I. INTRODUCTION 

The discovery of quasi-exactly solvable quantal problems was a remarkable contribution 
in last century[l]. This kind of problem could be solved by lie algebraic methods, which 
was pioneered by Turbiner[2]. The method was also be applied to an anharmonic oscillators 
by Chen et. al.[3]. Moreover a generalization of the basic lie algebraic construction was 
considered by Shifman[4]. You can also confer to a reviewed paper by Shifman[l]. However, 
the quasi-exactly solvable problem could also be calculated by an analytical way [5]. 

This paper is concentrated on a model that is a two-dimensional hydrogen with a linear 
potential in a homogeneous magnetic field. This model is very meaningful among the quasi- 
exactly problems. First, its potential is novel and complex, moreover when the coefficient 
of the linear potential becomes zero, its all result could go back to the simpler circumstance 
which have been investigated by Taut [5]. Second, the si (2) algebra structure hidden in the 
Hamiltonian is not apparent, however it will be revealed in this paper. Third, this model is 
also solved by an analytical method whose result coincide with the one by the si (2) method. 
Last but not least, a new sextic oscillator with a centrifugal potential is obtained from this 
model. Furthermore, the connection between them is constructed. 

This paper is organised as follows: in section two, the separation of variables in the 
corresponding schrodinger equation is displayed, so we get the radial part equation; in 
section three, we solve it by the lie-algebraic methods, i.e., sl(2); in section four, it is solved 
by an analytic way, therefore, the two results coincide with each other; in section five we 
build its connection with the sextic oscillator with a centrifugal barrier; in section six, a 
conclusion is drawed. 

II. THE SEPARATION OF VARIABLES 

Now, we consider a two-dimensional hydrogen with a linear potential in a homogeneous 
magnetic field, where the magnetic field B is perpendicular to the plan in which the electron 
is located. The schrodinger equation reads 

1 1 ^ Z 

[-{p + -Ay + kr]ilj = Ei/j , 

2 c r 

where c is the velocity of light and the vector potential in the symmetric gauge is .4 = ^Bxf. 
Under some operation, the above equation can be transformed to be 



1 1^2^ 1^ Z 

-\v'^i) + -(p ■ A)^ + -^- vi^ + ^^^V'l — ^ + ^fi) = E^ 

2 c c c^ r 



(1) 



In polar coordinates, letting r and 6' represent the radius and azimuthal angle respectively, 
we can get 



P= -He^^ + 



1 a 



A=\B X r = \Bre6 



A-p 



1 D d 



-t^B 



A^ = i5V2 



Substitution the above equations into the Eq.(l), we get 



1 d Z 

2 ov r 



(2) 



ld_^d_ I 1 a 

r ar dr r^ 89"^ 



and f^L = ^- Having represented the equation in a concrete 



where V^ 

coordinate, we will separate the variables in the next paragraph. 
Substituting 



Hr,0) 



e'"''R{r) 



(3) 



into the above equation, where m = 0, ±1, ±2, ..., we get the radial part equation, which is 

(4) 



1 d^ 11 d ulr'^ , Z m\^ ^^ 



The above differential equation will be solved by the method of sl(2) algebra in the next 



secion. 



III. SOLVING THE RADIAL EQUATION WITH THE SL(2) ALGEBRA 

First of all, let us recall the essential idea that the schrodinger equation is solved by sl(2). 
The schrodinger equation determining the stationary states of the system is 

Hi/j{r) = Eijj{r), 



According to [1], we can perform a quasi-gauge transformation 



iplr) = il){r)exp{— I A{r)dr). (5) 

Then the stationary schrodinger equation can be transformed to 

Furthermore H'^g" can be represented by the combination of the partial algebraization of 
the generators of sl(2), which is 

i7„c„ = J2 CabT'^T' + Y, CaT" + const, (6) 

where 



a,b-- 


=±,0 




a- 


=±,0 




r+ 


= 2jr 


— r 


2 ^ 

dr 




T^O 


T- = 


+ r 

d 

dr' 


d 
dr^ 



(7) 



It is very easy to check that the commutation relations for T^ and T^ are indeed those 
inherent to the si (2) algebra: 

[t^,t-] = 2t\ [T°,r+] = +r+, [t\t-] = ~t-. 

Finally, with the aid of the above commutation relations, the transformed Hamiltonian //"g" 
can be represented by a {2j + 1) x {2j + 1) matrix, whose basis are 

Then the standard algorithm of diagonalization of the finite matrix gives us {2j + 1) eigen- 
values and the corresponding eigenf unctions of the transformed Hamiltonian. 

Blindly following the program outlined above, we do not find the above model that we 
pay attention to in the section 2 can be solved by the methods of sl(2), in another word, 
the Hamiltonian H isn't reducible to the quadratic combination of the si (2) generators. 
However, we can consider a corresponding eigenvalue problem, /(r)iJ'0 = Ei(j, where / 



is an appropriately chosen function of the variable r, so that fH can be reducible to the 
standard form (6). The above explanation maybe a litter obscure, moreover this method 
will be showed in more transparent terms in the following paragraphs. 
Let us introduce an auxiliary Hamiltonian 



# = 7(- 



2d? 



l_d_ 
2r dr 



2^2 



UrT 



m 



2 2r2' 



+ ^kr — -jX, 



(8) 



where 7 and X are constant coefficients. Next, we calculate the "gauge'-transformed Hamil- 
tonian H"G", simply stemming form H after the substitution 



d d 
dr dr 



Then 



H,^G" = 7(- 



1 d^ 



2 dr"^ 2r dr 



Id d ^ 1 ^r, uir^ , 1 ^, A m 

+A—) — 7^^ +7—^^ V-fkr^ — 7A+7 h7 



Z 



dr' 



2r 



2^2 J, 



(9) 



+'^{muJL-X). 
(10) 



After a short analysis, we can get a conclusion that: when r — )■ cxd, the term 7—^ — V^kr will 

2 

become the dominant one; when r — > 0, the term 7^ will become the dominant counterpart. 
In order to guarantee the normalizability of the wave function, we must choose a proper form 
of A to cancel out the above terms. A short reflection shows that 



A = iir + 5+ -. 
r 

Substituting the above Eq.(ll) into Eq.(lO), one gets 



(11) 



+|7(m2 - z/2)^ + 7(+WLm ^ y. - iiv - \b'^ ~ X) 

(12) 

Moreover, from above analysis, one knows the coefficients of r^, r and ^ are equal to 

respectively, which are 

A; — /i(J = 



The above equations can be solved out: 



\i = ±Ul 



\ S 



(13) 



u = ±|m| 



Substituting Eq.(ll) and Eq.(13) into Eq.(14), one finds that only 

~ I k 
R{r) = R{r) exp( — OJ^r^ r + \m\ Inr) (14) 

2 lol 

can guarantee that the wave function converge at both sides of r — > and r — )■ oo(Note that 
we replace iplr) and iplr) by R{r) and R{r) respectively.), that is to say 

S=A. . (15) 

u = — |m| 
Substituting Eq.(15) into Eq.(12), one can gets 

Lh„^„r = E^r^ (16) 

7 7 

where 

+ [7;:^(l"^l + -2)-Z]+ 7K(1 + m + |m|) - 1(^)2 - X]. 
The reason that I divided both sides of Eq.(16) by 7 is for a convenient calculation. The reader 
maybe puzzled by the reason why we use r to multiple H"g"j because rH"G" can reducible to 
a quadratic combination of the generator of sl(2). Under the inverse transformation Eq.(9), 
the Eq.(16) can be transformed to be 

^HR{r) = ^^R{r), (17) 

7 ■yr 



where 



R=Rexp[ — ULr r+\m\mr). 

2 U>L 



Substituting Eq.(8) into Eq.(17), we can get 

^ 1 d^ Id ulr^ , IZ + E m\^ ^^ ^ 

^--2ir^-TrTr^^^ "- + """ - -,^- + 2;5)« = ^«- (18) 

By comparison with the Eq.(4), we observe that we have actually solved the original problem 
provided that the parameter 7 is chosen as follows 

^±1 = Z. (19) 

7 

And we represented the eigenvalue by another symbol X instead. 
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In order to solve Eq.(18) by sl(2), let's return to Eq.(16). Following what I said before, 
the -H"G" can be reducible to the combination of generators of sl(2), which is 



r ~ 



-H^'G" = CiT^'T- + C2T+ + CsT"" + C4T- + Co. 



7 
Substituting Eq.(7) into the above Eq.(20), we can get 

-H,,G" = C,r^ - C^r'^ + C,r^ + (Q - jC^)-^ + 2C2jr - jC+C,. 
7 dr^ dr dr ar 

By comparison with Eq.(16), the coefficients can be solved from equations 



(20) 



Ci = 


1 

2 


-C2 


= ujl 


Cs = 


_k_ 


c,- 


jCi- 



\m\ 



2C2J = ujl{1 +m+ \m\) - X 
-jC2 + Co = :^i^, + \m\)-f 

A further simplified calculation show that 



C2 



1 
"2 

-UJl 



^L 



C4 = -|(l+j + 2|m|) 

X = UL{2j + l + m+\m\)-U^f 



^o = ^(H + i + j)-f 



Substituting the above Eq.(21) and Eq.(20) into Eq.(16), one gets 



(21) 



[_1t°T- - ulT"- + —T' - ^(1 + J + \m\)r 



k n , 1 M ~ Z + E ~ , , 

— {\m\^-^^3)\R = R (22) 

Wl 2 7 



By use of the constraint which is Eq.(19), the above Eq.(22) can be transformed to 

[_1tOt- _ ^^T+ + Aro - 1(1 + J + \m\)T- + A(|^| + 1 + jy^R = ZR (23) 

Z UJl Z UJl z 



i)j=o 



Hence, we know 

Substituting Eq.(24) into Eq.(14), one get 

1 



R = ttQ. 



k 



i? = aoexp( — ult r+|m|lnr). 

2 ul 

Substituting Eq.(7) and Eq.(24) into Eq.(23), we can get 

-{\m\ + h = Z, 
ul 2 



(24) 



(25) 



(26) 



which is precisely the constraint of eigenvalues. Substituting the above Eq.(26) into the fifth 
Equation of Eq.(21), we can get the eigenvalue 

1 Z 

Xo = uiil + m + \m\) - -{- ^2 



2'|m| + |' 



2)J = I 

It is very easy to calculate that 

1 




jiO 



1/10 
2^0-1 



T' 





1 



(27) 



Substitution the above Eq.(27) into Eq.(23), we can get 



^(|m| + |) -UL 

-(\m\ + 



— imi 



1> i^L 



ai 
ao 



«i 



ao 



(28) 



Such a equation has a nontrivial solution if the determinant of the corresponding matrix 
vanishes, i.e. 

-(H + 1) ^(|m| + i)-Z 
The equation can be further simplified to be 



0. 



ujIZ-" - 2kujLi\m\ + 1)Z + k^i\m\ + ^)(|m| + 1) - a;i(|m| + 1)] = 0. 



(29) 



Form it, we can draw a conclusion that the parameters Ul and k can't be chosen arbitrarily, 
they must satisfy the above Eq.(29). The above Eq.(29) can be regarded as a quadratic 
Equation about Z. As a result, the solution of it can be expressed as 

k{\m\ + 1) ± ^/M 



ojl 



(30) 



where M = ^ + Ui{-^ + \m\). In the following paragraph, we will calculate the different 
eigenf unctions corresponding to different eigenvalues. 

g\2 = k{\m\+l)+VM 



^L 



Substituting the above equation into Eq.(28), we get 



ao 



1^ + ^VM 

1 



«o(l - 



lA + ^^M 

i + \m\ 



r . 



(31) 



Substituting the above Eq.(31) into Eq.(14) 



J 2 



2 Wi ijJL ^ 



1 k 

-I r)exp( — (jOhf"^ r+|?72|lnr). 

k + \m\ 2 ujl 



fc(|m|+l)-VM 



b)Z 

Following the above procedure, we get the radial part of the wave function, i.e. 



l^L 



R{rr: 

J~2 



ao(l 



i-^-^VM 1 k 

-rj exp( — UlT r+|m|inrj. 



1 



\m\ 



2 , i-.^i - ^L 

With the constraint of Eq.(29), the eigenvalue of energy becomes 

1 k 

Xi =ujL{2 + m+\m\)--{ — f 

3) Arbitrary j 

According to the choice of the phase, we know the following formulation, i.e. 



(32) 



(33) 



T \jm) = Y (j =F m)(j ± m + l)|j, mil) T^ljm) = m\jm). 
The j± and Jq can be represented in the following {2j + 1) x (2j + 1) matrix. 



T' 



10 VI- 2j 






\0 


v^ 











\ 

1 
J 


(2j- 






-1) 


■ V(2j-l)-2 
V2j • 




(34) 



T- 



( 




V 








v/2~(27 



o\ 







V(2j - 1) ■ 2 



(35) 



/ 



T->U 



v/2jTT oy 



J-1 



J -2 



-J + 1 



V 



(36) 



-^7 



Substitution the above Eq.(34), Eq.(35) and Eq.(36) into Eq.(23), so the problem is deduced 
to an eigenvalue and eigenfunction problem. In order to get the nontrivial eigenfunction, 
the determinant of the corresponding matrix must vanish, which will give a onstraint about 
bJi^ and A;. Similarly, we will get the radial part of the wave function. Furthermore, the fifth 
equation of Eq.(21) tell us the eigenvalue of energy. 



X,=a;i(2j + 1 



va 



\m\ 






To sum up, in this section, we have solved Eq.(4) by method of sl(2). Moreover, we will 
solved it again by another method, namely, an analytical methods. 



IV. SOLVING THE SAME EQUATION WITH ANALYTICAL METHOD 



In order to solve Eq.(4), using the following transformation, i.e. 



we can get 



l(fu 1 



1, 1 



[m 



:U 



R — —j=, 



2 2 7-^ 

—u + Utttiu + kru H — UJrr u = Eu. 
r 2 ^ 



(37) 



(38) 



2dr2 2' 4V2 

Due to the condition of square-integrability of the wave function, it is useful to consider of 
the limits r — )■ and r — )■ oo of the above Eq.(38). For r — ?■ 0, the term containing \ is 



10 



dominant and we get the equation 

1 (fu 1 , o 1 X 1 / s 

Because the above Eq.(39) contains a normal singularity, so we substitute a power series 
'" = ^"^ Si^o ^«^* ^^^'^ ^^^ equation and one gets 



i=0 j=0 



^s+i-2 _ Q^ 



Due to the uniqueness of the series, the coefficient of it is equal to zero, we can get 

-(s + i){s + i- l)bi + (m^ - -)bi = 0. 

By use of the limits r — )> 0, only the lowest-order term become dominant, so we omits the 
higher order terms and get 

-s{s - l)bo + K - 4)^0 = 0. 

In order to get the nontrivial solution, we must set bo ^ 0, so one gets s = m + | or 
s = —m + i. However, if m become negative, the solution diverges when r — )> 0. Hence, we 

get 

, , 1 

s = \m\ H — . 

So the solution of Eq.(39) is 

n = rH+5. (40) 

Having solved the solution r — !■ 0, we consider about another situation r — )■ 00. Under this 
circumstance, the terms containing r^ and r in Eq.(38) become dominant. So the equation 

becomes 

1 cPu 1 2 2 7 

-— H — ujrr u + kru = U. 

2dr^ 2 ^ 

We can verify the solution of the above equation is 

M = exp(— -w^r^ r). (41) 

2 ujl 



Considering both Eq.(40) and Eq.(41), we can write the solution of Eq.(38)as 

^^^IH+lg-^^^^'-^V. (42) 

Substituting the above Eq.(42)into (38), we get 



-|r$ + u.r^'j + ^r J - (|m| + 1) J + {\m\u, + u, + mu, - l| - E)r^ 

11 



(43) 



The above differential equation has two nonessential singularities at r = 0, oo. It's index 
equation takes this form 

s(s + 2|m|) =0. 

It's solutions are and — 2|?7i| respectively. However, if s = — 2|r?2|, the solution diverges, in 
another words, it violate the demand of square-integrability of the wave function. Hence, s 
must be 0. Therefore, we can substitute the power series 



^ = J2<^nr'' (44) 



n=0 

into Eq.(43)and find the recurrence relation 

Let us discuss in detail in the following paragraphs. 
l)n = 

ai = I ,""i ao. (45) 

We put forward a hypothesis that the ground state exists, so we get ao 7^ 0, ai = 0, 02 = 0, 

■ ■ ■ . However, conferring to Eq.(46), if Oq 7^ 0, a2 may be not equal to 0, which is contradict 
the hypothesis. So we can draw a conclusion that oq = 0, that is to say, there isn't a ground 
state. 

2)n = 1 

\m\ujL + ujl + muL - l^ - E nm\ + ^)^- Z 

02 = 7, — ; ^ ao H 7, — ; — ^ CLi- 46) 

2(|m| + l) 2(|m| + l) ^ ' 

If the first excited state exists, ao and ai can't be equal to simultaneously, 02 = 0, as = 0, 

■ • ■ . Conferring to Eq.(49), ai will have a effect on 03, so ai = 0, that is to say, conferring 

to Eq.(45), 

Z={\m\ + \) — . 
2 ujl 

So the radial part of wave function is 

Rn=i = rl-le-^^^^'-^'-ao, (47) 



12 



El = ujl{1 + \m\ + m) - - — . (48) 

Z Ut 



which is gained by use of Eq.(37) and Eq.(42).By comparison with Eq.(25), it is found easily 
that the two results coincide with each other. From above analysis, we can draw a conclusion 
that ao 7^ 0, ai = 0, a2 = 0, ■ ■ ■ , conferring to Eq.(46), so we get 

A 

3)n = 2 

3(|m| + |) 3(|m| + |) ^ ^ 

Following the above analysis, if the second excited state exists, we can get 03 = 0, 04 = 
and so on. It is necessary to demand the first term and the second term are equal to 
respectively, conferring to Eq.(49). Moreover if a2 7^ 0, it may render 04 7^ 0, which is 
contradict our hypothesis. So we must let 02 = 0, in another word, the second term become 
zero already. Now let us concentrate on the first term. If ai = 0, this case will be the same 
with the first excited state. So we can get 

1 k"^ 
E2 = {2 + \m\+m)LOL- -^. (50) 

2u^ 

Substituting Eq.(45) and Eq.(50) into Eq.(46), we can get 

_ [(l"^l + 1)4 - ZEH + 1)4 -Z]- (|m| + IP, 
2(|m| + l)(|m| + |) 

From above analysis, we know 02 = 0, that is to say, 

[iH + h- - Z][{\m\ + h- -Z]- {\m\ + hu, = 0, 

which coincide with Eq.(29). Therefore, regarding Z as unknown root, we can certainly get 
the same solution with Eq.(29), i.e., Eq.(30). Substituting Eq.(42) and Eq.(44) into Eq.(37), 
we can also get the radial part of the wave function, 

Rn=2 = r' 'e 2 ^ -i ■ — ■ — k r + 1 ao- (51) 

\m\ + i 

Substituting Eq.(30) into the above Eq.(51), we get 

\m\ + 2 
which is coincide with Eq.(32) and Eq.(33). 

13 



4:)n>2 



cuLin + \m\ + m) - ^fjT - E ^\m\ + ^ + n)^- Z 

{n + l){\m\ + ^) """-' + (n + l){\m\+ ^ 



O-n+l — -, ~TT7i 1 , ii„N Cln-l H } "TTTl 1 , i+„x '^n- (52) 



En = u}Lin+\m\+m) - -^, (53) 

z ujt 



Following the similar analysis above, also conferring to [5], it is a sufficient condition to 
guarantee normalizability of the eigenfunctions, i.e., 

A 

and 

dn = F{\m\,n, E,ujl, k, Z)ao = 0. 

For ao 7^ 0, we can deduce 

F{\m\,n,E,UL,k,Z) = 0. (54) 

Substituting Eq.(53)into above Eq.(54), we obtain 

1 k"^ 
F{\m\,n,ULin + \m\ + m) - -^,UL,k,Z) = 0. (55) 

From it, we obtain a constraint for ul and k, in another word, ul and k can't be chosen 
arbitrarily. Furthermore, from Eq.(45), Eq.(46) and Eq.(52), we can get the polynomials of 
ip, i.e. Eq.(44). By use of Eq.(37) and Eq.(42), we can get the corresponding Rn- 

In conclusion, we have solved the Eq. (4) by an analytical method. In the following section, 
we will build the connection between the radial part equation and the sextic oscillator with 
a centrifugal potential. 

V. KS TRANSFORMATION 

Via KS transformation, i.e., 

2 

r = p 

(56) 

9 = 2^ 
the Eq.(2) becomes 

2 [- V (a^ + -pd-p^7W^^- ^"^¥^ + ^^'^^ - 7^ + ''^ = ^^- 

Substituting 

^(p,^) = ^e^'^^x(p) (57) 

V 27r 
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into the above equation, where m = 0, ±1, ±2, ..., we get 

1 d 1 d Tfi 

2 dp^ prfp p^ 

In order to cancel the item which contains -, substituting 

X = P"k (58) 

into the above equation, one obtains 

which is the so-called sextic oscillator with a centrifugal barrier [6]. By comparison between 

Eq.(3) and Eq.(57), we can get 

in 
m = -, (60) 

and 

x{p) = R{r). (61) 

Substituting Eq.(58) into the above Eq.(61), we obtain 

C = v^i?(r), (62) 

which is precisely the wave function of the sextic oscillator with a centrifugal barrier. More- 
over, using Eq.(56) and Eq.(60), we can replace the parameters in Eq.(62) by the coun- 
terparts used in Eq.(59). In section two and three, we have calculate Z by two different 
methods. Therefore, it is very easy to get the eigenvalue 4Z of Eq.(59). 
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